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Abstract 

The toxins associated with infectious diseases are potential targets for inhibitors which 
have the potential for prophylactic or therapeutic use. Many antibodies have been gen- 
erated for this purpose, and the objective of this study was to develop a simple mathe- 
matical model that may be used to evaluate the potential protective effect of antibodies. 
This model was used to evaluate the contributions of antibody affinity and concentra- 
tion to reducing antibody-receptor complex formation and internalization. The model 
also enables prediction of the antibody kinetic constants and concentration required to 
provide a specified degree of protection. We hope that this model, once validated exper- 
imentally, will be a useful tool for in vitro selection of potentially protective antibodies 
for progression to in vivo evaluation. 
Keywords: toxins; antibodies; kinetic model 

1. Introduction 

Passive immunization using antibodies has been used successfully for treatment and 
prophylaxis of infectious disease in humans [l| and there is increasing interest in the use 
of antibodies for treatment of infectious diseases that may be used as terrorist weapons, 
but for which the risk is not sufficiently high to justify preventive vaccination of a large 
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civilian population [2] . Toxins are an important potential target for designing therapies 
against these threats and a broad range of approaches has been taken to develop inhibitors 
that may be of prophylactic or therapeutic use jjj. 

Antibody engineering techniques allow affinity maturation of antibodies and these 
techniques are being exploited to produce inhibitors for a number of toxins [jj]. The 
emphasis of this approach is on producing reagents with high affinity, based on the 
proposition that higher affinity will provide better protection. 

However affinity, by itself, is a poor predictor of protective or therapeutic potential. 
Antibodies with high in vitro affinity for toxins do not automatically confer protection 
in vivo [fjj], [3] and may exacerbate the toxicity [gj]. The effects of using multiple 
antibodies with high affinities may be additive 1(3] or synergistic or without effect 



In addition, epitope specificity 



llj . antibody titre [l2| - [l6jj and dissociation rate 



have been correlated with protection. 

Toxins are produced by a number of plants, animals and microorganisms. Toxins may 
act at the cell surface and either damage the cytoplasmic membrane or bind to a receptor 

n 

and act via transmembrane signalling subsequent to that binding 11811 - Alternatively, 
toxins may cross the cell membrane and act on intracellular targets [la ]. For example: 
anthrax lethal toxin, ricin and cholera toxin bind to a cell surface receptor and make use 
of cellular membrane trafficking to enter the cell [l^ . [20I ]. 

The objective of this study is to develop a simple mathematical model that may be 
used to predict the optimum antibody parameters (kinetic constants and concentration) 
needed to inhibit the binding of the toxin to its receptor. These predictions may be 
used to select candidate antibodies for progression to in vivo evaluation and to assess the 
potential value of affinity enhancement. 

This paper is an extension to our previous work 21]. In the model presented below 



we explicitly take into account the process of toxin internalization and diffusive fluxes 
around the cell. 



2. Model 

The kinetic model describing the interactions of toxins with cell receptors can be 

formulated based on the well-known analytical framework for ligand-receptor binding. 
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Figure 1: Schematic representation of the model for receptor-toxin-antibody interaction. 

The models of this process have been studied for many years and a vast amount of 



22i- m and refs therein) 



literature has accumulated on this subject (see 

When a toxin diffuses in the extracellular environment and binds to the cell surface 
receptors, the toxin concentration will vary both in space and time. Any rigorous descrip- 
tion of this process would entail a system of Partial Differential Equations (PDE), which 
couples extracellular diffusion with reaction kinetics of the cell surface. The resulting 
system of PDE is nonlinear and too complex to be treated analytically. This complexity 
makes unfeasible any comprehensive study of parameter optimization. From another 



22j - |26| and 



perspective, it is well known that under some rather broad conditions (see 
refs) the reaction-diffusion system of the ligand-receptor binding can be well approxi- 
mated by a system of Ordinary Differential Equations in which the spatial variability of 
the process is simulated by different concentrations of species in initially predefined spa- 
tial domains (called compartments). Although, this compartment model is significantly 



simpler than the initial reaction-diffusion system, it still allows a consistent description 
of diffusion fluxes [23 1, Q, 2fl- 



In the current paper we use the compartment-model 
approach for our analytical study and numerical simulations. 
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To begin, we consider the following simple model. The toxin, T, binds reversibly to 
cell surface receptors, R, with a forward rate k\ and a reverse rate to form the toxin- 
receptor complex Cr which is then slowly internalized at a rate k% . The neutralizing 
antibody binds competitively to the toxin with on and off rates of k 2 and fc_2 respectively. 
The antibody-toxin complex, Ca, remains in the extracellular space (see FigfT]). 

We can easily write an equation for the toxin- receptor binding (viz, without antibody 
being present). For a spherical cell of radius a with the toxin binding to its surface (22I ]- 

H 

rfC 

^ = k e f RT + k e r C R) (1) 

where Cr is the concentration of the bound receptors (toxin- receptor complexes), R is 
the concentration of receptors, T is the bulk toxin concentration (i.e. far from the cell 
surface) and is assumed to be spatially uniform. The effective forward and reverse rate 



coefficients are defined by 



22 



2G] 



k) = 7 fci, K=lk-i, (2) 
where ki,k-\ are intrinsic reaction rates, kjj = AitaD is the diffusion rate, D is the 
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diffusivity of toxin in the extracellular space, 7 = 1/(1 + Rki/kjj) < 1 [23j.|24|. 

The bulk concentration of toxin T is mainly driven by the binding to antibody. There- 
fore, in this case we can write an equation system similar to (|Tj), but without any "dif- 
fusive" modification of the intrinsic rate constants: 

= k 2 AT + k. 2 C A , (3) 

where Ca is the concentration of toxin-antibody complexes, A is the concentration of 
antibody. 

The process of toxin internalization is phenomenologically introduced into our model 
by the equation 

^ = hC R , (4) 

where Ti is the concentration of internalized toxin. The corresponding term should be 
included in ([1}, so we arrive at modified expression for k^ 

fc^=7fc_i-/c 3 . (5) 
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The system (JTJ) , ([3]), (j4]) should be supplemented with three conservation laws for 
concentrations of R, T and A: 



R = R + C R , (6) 
A Q = A + C A , (7) 
T Q = T + C T + C A + T l , (8) 

where Rq, To and Aq are the initial concentrations. 

Eqs (HJ), ([3]), (O - © form a framework for our analysis. This is a system of 
nonlinear ODE (because of conservation laws ([6]) - ([8]) and because of effective rates 
k*j , are functions of the receptor concentration) . It can be easily solved numerically 
and also allows some analytical progress (see below). If parameter 7CI (and this is the 
case in many practical situations), then this model can be reduced to the "well- mixed" 
kinetic model with constant kinetic rates 21 1 . 



It is worth emphasizing that the aim our analytical framework is to develop a simple, 
but scientifically rigorous model that may be used to predict the optimum antibody ki- 
netic properties and concentration required to achieve a desired protective effect, rather 
than develop a detailed, biologically accurate model that captures all the detail of the 
toxin internalization process. Therefore, the model does not take into account the phar- 



macokinetics of the toxin-antibody complex |9| or receptor internalization and recycling 
2^ |. [2^. &3 is a lumped constant that describes all the processes that result in the ap- 



3l| and 



pearance of the free toxin in the intracellular space [3(3 ] . Wiley and Cunningham 
Shankaran et al. 32| have also developed mathematical models of this type of process. 

We are particularly interested in the behaviour of the model under conditions most 
likely to reflect the real biological situation i.e. toxin concentration much lower than the 
concentration of receptors (Tq/Ro <1). 

Testing of the model was carried out using COP AS I (software application for simulation 



331 ]) and the kinetic parameters 



and analysis of biochemical networks and their dynamics 
for the binding of ricin to its receptor and its internalization [34( and competition by the 
monoclonal antibody 2B11 6]. The kinetic parameters used are shown in Table 1. The 

m 

value of &3 used is that determined by Sandvig et al. [19( to be the rate of irreversible 
binding of ricin to HeLa cells. For simplicity, the simulation was carried out using all 
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Table 1: Kinetic constants used in numerical simulations (the binding of ricin to its receptor and the 
monoclonal antibody 2B11). 

Reaction Value 

fci 1.3 • 10 5 

fc_i 1.4 ■ 1CT 2 

fc 2 1.25 • 10 5 

fc_ 2 5.2 ■ 10~ 4 

fc_ 3 3.3 -10- 5 



reactions taking place in the same compartment. 

To illustrate the model we used toxin and receptor concentrations based on cell culture 
studies carried out in our laboratory. These typically use a cell concentration of 1-10 4 cells 

)er 100 /il experiment and a ricin concentration of 10 pM. Assuming 3 • 10 7 receptors/cell 

34^ . the receptor concentration is approximately 5nM. 



3. Analytical Results 

3.1. Cell- Surface Binding 

Initially we derive some analytical results for toxins that act at the cell surface and 
are not internalized, i.e. we set — in Eq. Q. At equilibrium d/dt = and from if!]). 
([3]) we can write 

C R = RT/Ki, C A = AT/K 2l (9) 

where K\ = ki/k—i, K 2 = k 2 /k-2 are the association constants for the toxin binding to 
the receptor and antibody respectively. It is worth noting that the parameter 7 (diffusive 
correction of the intrinsic reaction rates) disappears from Eqs.©, so in this case the 
analytical results are identical to ones derived using the "well-mixed" approximation 
21|. 

In order to simplify notations we denote by z and y the equilibrium concentrations 
of the toxin-receptor and toxin-antibody complexes, i.e. 

Z=[Cn]e q , V=[C A ]e q - (10) 
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From Eqs.(|9]) and conservation laws ([6])- ([8} the following closed equation for z can be 
derived: 



(R -z){T -z-y)-K 1 z = 0, (11) 

cz 

Ro - z{l - e) 

where e = K\jKi- 

Eq. ljlip can be written in a more conventional form of a cubic equation: 

a3z 3 + a2z 2 + a\z + ao = 0, (13) 

where 

a 3 = e - 1, 
a 2 = (l-e)C + eA + i?o,r 
oi = -i?o(Co + A + (l-e)T ), 
»o = ^o-Rq, 

and Co = Ro + K\ . 

It is well-known that Ea. dTBl has a closed- form analytical solution (Cardano's formula 
:.35|), which in our case provides a consistent way to derive exact solutions for the pro- 
posed model. Unfortunately these solutions still involve rather cumbersome expressions, 
which require further simplifications in order to be used in practical situations. Below 
we present another approach that explicitly employs the smallness of ratio Tq/Rq <C 1 
and leads to a simple analytical expression for the protective properties of the antibody. 

We observe that in the absence of antibody (i.e. Ao = 0), Eq. (fTTj) is an elementary 
quadratic equation that has two roots. If we impose the obvious constraint z —¥ as 
To — > then there is only one solution, which we designate as zq: 

Co |, ( A 4i? 7o x 





*-f '-l 1 — arl (14) 



Under the condition Tq/Ro <C 1 , this can be simplified to 



zo~^, C = Ro + K 1 . (15) 
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Let us now evaluate the effect of adding an antibody. From a mathematical point 
of view this effect (i.e. change of z under condition A$ > 0) is captured entirely by the 
term y in so our aim is to provide a reasonable analytical estimation of this term. 

From Eq. (I12j) and based on our initial assumption of low toxin concentration (To /Rq -C 
1) we can deduce the following simple estimate y ~ ezAo/Rg. This then leads to a mod- 
ified form of Eq. (fTTj) 

(Ro-z)(T -z)-K*z = 0, (16) 

where 

K* = Ki + eA . (17) 

We can see that this is the same form as the equation for z when Aq = 0, but now 
with K\ replaced with K*. This also implies that the analytical solution (|15|) is still 
valid, but only with the substitution K\ = K*. 

In order to characterize the effect of an antibody on the binding of a toxin to its 
receptor, we introduce the non-dimensional parameter the relative reduction in Cr 
due to the introduction of an antibody 

z(A Q > 0) 

z(A = 0) { } 

The analytical results presented above enable us easily to derive a simple formula for the 
antibody efficiency parameter By using (fTU]) . [T5l (IT71) . (|T%|) we can readily deduce 

*=-J— , e = K 1 /K 2 , X = A /C . (19) 

1 + 6A 

This expression is the main result of the current paper and will be validated with nu- 
merical simulations. 

To conclude this section let us briefly discuss some additional constraints for the 
parameters of our model in order for the expression (1191) to be valid. As mentioned 
above the condition of low toxin concentration is always assumed in our study. Another 
simple condition can be derived from the constraint Cr + Ca < Tq and by using Eq. (1151) : 

£(l + e£)«4£<l, (20) 
Oo Oo o 

since Rq/Co < 1. This condition could be always checked retrospectively and always 
hold in our numerical simulations. 



3.2. Toxin Internalization 

For toxins that are internalized, the effect of antibodies that prevent receptor binding 
is to reduce the effective rate of internalization. To examine and evaluate this effect we 
need to analyze the full system ([1]), ([3]), ©. 

In order to characterize the effect of antibody concentration on the rate of toxin 
internalization we introduce a new parameter: 

which is a function of time (i.e. G = G(t)). 

Our aim is to deduce function G based on the kinetic model (JXJ) , © , © ■ It is 
evident that G < 1 for t > and G — > 1 as t — ► oo (since in that case all toxin will be 
internalized) . 

For the toxins of interest, while the receptor binding is rapid (time sale ~ l/(fciCb)) 
22 1, [23I ]. the subsequent internalization is much slower (time scale ~ 1/^3 ^> l/(fciCo)). 



This coupling of slow and fast processes in our system allows us to develop a simplified 
model of toxin internalization usin g th e the well-known framework of Quasi-Steady-State 



Approximation (QSSA), see 22|- [27] and refs therein. 

When applied to our system QSSA elucidates the toxin internalization as a two- 
stage process. After the initial rapid binding of the toxin to the receptor we can simply 
set dCn/dt = in ©. The further slow evolution of T(t) (viz, quasi-steady state) is 
completely determined by the conservation law © and Eq.© and spans a time scale of 
the order of the internalization time (~ 1/ks). In addition, for solving © at the initial 
stage of internalization, we can assume that T <C To and write 

Ti(t) = k 3 z t, i«l/fe 3 , (22) 

where zq is given by expressions (Q© and ([T5|) . The evolution of Ti(t) for the late stage 
of internalization can be readily derived from ©, © - © by assuming [To — T(i)] <C To: 

T,(t) = T [l - cxp(-M)], t > 1/ks, (23) 

so Ti(t) exponentially approaches its saturation limit. A simulation of this process is 

shown in Figj5] and the slow linear increase of T at the initial stage is clearly visible. 
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Now, consider the case where Aq > 0. According to (1221) the main effect of the 
introduction of an antibody is to reduce the value of zq, as described in the previous 
section. Then, based on (|2"Tj) . (|2"2"j) and (TT51) we can conclude that, during the quasi- 
equilibrium stage, the following approximation holds 

where if? is given by expression (|T5|) . 

The overall effect of introducing an antibody can be best describes in terms of the 
internalization half-time, r,. Without antibody the later can be estimated from (|2U)) and 
condition T^n) = T /2. Thus from Ea.([2"2"]) we yield 

~ — = ^° (2^) 
1 ~ 2fc 3 z 2fc 3 i? ' 1 j 

For the internalization time with the presence of antibody we can apply reduced value 

of zq and write the following simple formula 

5 * < 26 > 
where t 4 ° is the internalization time in the absence of antibody (Aq = 0). 

Eqs. ([2"5|) and (121)1) have a clear interpretation. As described in the previous section, 
the introduction of an antibody results in a decrease, at t <C Tj, in the equilibrium value 
of Cr (i.e. in zq)- This can be related, in accordance with Ea. (|22|) and Ea. (|25|) . to a 
corresponding decrease in the concentration of internalized toxin Tj and a consequent 
increase in the toxin internalization time (since it takes longer to achieve a give level of 
Tj). Since changes in z$ can be described comprehensively by means of the parameter 
it still remains the only parameter needed to characterize the influence of an antibody 
on the concentration of internalized toxin (IM|) . (|2l)|) . 

It is evident that the two main effects described above (reduction of the concentration 
of internalized toxin at a given time, and increase in the time required for the internalized 
toxin to reach a given concentration) are not independent of each other. The linear 
relationships ((M]), (121)1) allow us to establish a general identity that relates these two 
effects for any time t. 

Let us assume that for Ao — 0, tq in is the time taken for the internalized toxin to 

reach a concentration T? (i.e r° = Tf/{k^Zo), see d2U0- The effect of introducing an 
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antibody is to reduce the internalized toxin concentration to a value Ti < T® . Then 
from (j2"4")h (|2"r?)) we can derive the following identity: 

T z r t = T°r°, (27) 

where r, is the time required for the internalized toxin to reach TP when Aq > 0. The 
identity (|27l) has no explicit dependency on antibody kinetic parameters or concentration 
and provides an easy way to calculate any of the parameters (X^X^ , Tj,rP) if the other 
three are known. 

4. Numerical Results and Discussion 

We have derived an analytical expression for the parameter "J, the relative ability of 
an antibody to reduce the binding of a toxin to its receptor (|19p . Our derivation is based 
on the following assumptions: 

1. Toxin concentration is much lower than the receptor concentration 

2. For toxins that are internalized, the internalization rate is much slower than 
establishment of the receptor-toxin binding equilibrium. 

Applying these assumptions we found that parameter 'J is independent of the toxin 
concentration (see (fT9|) ). i.e. it is determined by the ratio of antibody to receptor con- 
centration and not by the ratio of antibody to toxin concentration as is commonly used. 
For the low toxin/receptor ratios likely to occur in biological situations, the condition 
(|20j) can be met by large range of antibody kinetic parameters. From this point of view 
Ea. (fl9l) should be valid for most practical applications. 

The implications of our analytical findings are illustrated by simulation of the com- 
plete kinetic model (Eqs. ((H,©, dU, ©) using the kinetic constants for ricin and 
the anti-ricin antibody 2B11 (Table 1). Figj2]is a simulation of the effect of the presence 
of an antibody on the binding of the toxin to its receptor (formation of Cr). The anti- 
body concentration is expressed as the dimcnsionlcss parameter A = Aq/C'q. In this case, 
since Rq and Tq <C K±, the parameter Cq = Ro + K\ is dominated by K\ (1.08 ■ 10 -7 ). 

Fig[3] shows the effect of increasing antibody concentration on ^ . There is a good 
agreement between the values of determined from (JTUJ) and from (jT5J) using the equi- 
librium values of Cr determined from simulation of the complete kinetic model (Fig. [3]) . 
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t (sec) 

Figure 2: Simulated effect of antibody concentration on formation of toxin-receptor complexes Cr. 
Parameter A = Aq /Co, Co = Ro + K\ . The binding curves were created using the simulation package 
C0PASI and the kinetic constants in Table 1. R = 5nM; To = lOpAf, Co = 1.15 ■ 10~ 7 . 
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Figure 3: Effect of antibody concentration on protection factor. Parameter ^ (1181) was determined from 
Eq. lll9l l (solid lines) and by using simulated values of Cb, from Fig[2]at 2500 sec (A), e = 25.9. 
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Figure 4: Protection factor dl8l as a function of parameter e = K1/K2 and A = Aq/Co (Eq J19l l): 
A = 0.01 (o); 0.025 (A); 0.05 (□); 0.1 (y); 0.25 (0). The range of values for A and e below dashed line 
corresponds to 80% protection. 




Figure 5: Different time scales for formation of receptor-toxin complex Cr (□) and associated toxin 
internalisation Xi (solid lines). Results of C0PASI simulation with kinetic constants from Tabic 1. A = 
A /Co, Ro = 5nM; T = lOpM, Co = 1.15 ■ 10~ 7 , e = 25.9. 
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3t 

Figure 6: Comparison of parameters and V. T (A) was determined using values of T; and T° at t = 10 4 
sec from toxin internalization time courses simulated using C0PASI and the kinetic constants in Table 
1. Parameter * (solid line) was determined from Eq. {l9)l . Rq = 5nM; T = lOpM, Co = 1.15 ■ 1CT 7 , 
e = 25.9. 

For instance, the results predict that, for this toxin and antibody combination, the addi- 
tional protection provided by increasing the antibody concentration diminishes rapidly 
when A exceeds 0.1. 

Fig0]shows the relationship (TIT)|) between antibody concentration and the toxin/antibody 
and the ratio of toxin/receptor dissociation constants (e). This plot is valid for all com- 
binations of toxin, receptor and antibody consistent with the assumptions used to derive 
(j!9[) . principally Tq <^ Rq. The antibody kinetic parameters and concentration required 
to provide a specified degree of protection may be determined from this plot. For exam- 
ple, any combination of e and A falling below the dashed line will reduce either Cr or Tj 
by 80 %. 

This, in turn, enables important judgements to be made about antibody selection. 

For example, if an antibody concentration of 0.25Co (A = 0.25) is achievable, then an 

antibody with an e value of 50 will provide good protection (93% reduction in Cr or Tj). 

If an antibody concentration of only 0.05C*o (A = 0.05) is achievable, then an e value of 

250 is required to achieve the same level of protection. The structure of (flT)|) is such that, 

a given increase in protection or T) may be achieved by either an x-fold increase in e 
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Figure 7: Establishment of the quasi-equilibrium state in the presence of antibody. Cr formation (A) 
was simulated using C0PASI and the kinetic constants in Table 1. T (□) was determined using Eg . (I24D 
and values T and TP at t = 10 4 sec using simulated toxin internalization time courses. Rq = 5nM; 
T = lOpM, C = 1.15 ■ 1CT 7 , A = 0.05. 



or an x-fold increase in A. 

The effect of antibody on toxin internalization is simulated in FigfSJ Rapid equili- 
bration of receptor and toxin is followed by slow accumulation of toxin within the cell. 
Equation Cj24]) predicts that VP is the only parameter needed to characterize the influence 
of an antibody on toxin internalization. FigJB] compares T calculated using (pM|) , (fT5|) 
with r determined using values of Ti and at t = 10 4 sec from this simulated data and 
shows good agreement between the two values under the condition To <C Rq, although 
the value of T is slightly greater than 'J. The plot predicts the degree of protection 
provided by a given concentration of antibody and enables assessment of the value of 
increasing antibody concentration beyond a certain value. For example, to enhance the 
reduction of T from 90% to 95% requires doubling of Aq. 

The expression for 'J, (|19p . assumes a quasi-equilibrium state in the system. In 

practice, this state may take significant time to achieve. Fig. [7] shows a simulation of 

the time taken by the ricin/receptor/2Bll system to reach the quasi-equilibrium state 

for A = 0.05. The value of T determined from the toxin internalization profiles (FigJT]) 

parallels this process, i.e. experimental validation of T must allow sufficient time to 
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Figure 8: Relationship between toxin internalization time tj and protection factor ^ H18II . Solid line is 
formula (126 I t and (□) is simulation with CQPASI. r; was determined as the time to internalize 5 • 10~ 14 M 
of ricin. All other parameters are the same as in Fig[7] 

elapse for the quasi-equilibrium state to be established. 

The relationship between the internalization time Tj and ^ described in Ea.(|26p is 
shown in Fig|8l VP was determined from simulated toxin internalization time courses 
(FigJS]) as the time to internalize 5 • 10 _14 M ricin. The slope of the fitted line is 1.07, 
close to the predicted value of 1.0. 

In summary, the protection provided by an antibody against toxins that act either 
at the cell surface or after binding to the cell surface followed by internalization may be 
predicted from a simple kinetic model. Protection parameter ^ is a simple function of 
antibody, receptor and toxin concentrations and the kinetic parameters governing the 
binding of the toxin to the receptor and antibody: 

* = l + ^/tfaXAo/Co)' (28) 
The calculated value of $ matches closely the degree of protection determined from 
numerical simulation of the binding and internalization reactions and provides a con- 
venient method for predicting the optimum antibody parameters (concentration and 
dissociation constant) needed to provide effective treatment or prophylaxis for toxins. 
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